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ABSTRACT 


A  means  is  described  for  computing  the  peak.  Static 
and  peak  dynamic  overpressures  of  the  shock  wave 
generated  by  an  artillery  piece  with  muzzle  brake  as 
a  function  of  position  within  the  crew  area.  Computed 
results  are  compared  with  recent  experimental  results 
and  indicate  a  favorable  agreement. 
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LIST  OF  SYMBOLS 


The  reader  is  referred  to  RIA  report  Number  62-1794 
for  a  complete  list  of  symbols  concerned  with  gun  interior 
ballistics  and  muzzle  brake  parameters.  A  partial  list 
is  repeated  here  for  reference. 


INTERIOR  BALLISTICS  SYMBOLS 


Symbol  Meaning 


Dimension 


A 


^cnct 


Mq 


area  of  bore  plus  cross  sectional  area 
of  grooves 

area  of  ammunition  band  in  contact  with 
tube  after  shear  deformation 

total  internal  volume  of  gun 

mass  of  projectile 

mass  of  charge 
mass  of  igniter 

total  propellant  mass 


in2 


lbm 

lbm 

lbm 


Mr  mass  of  recoiling  parts 

T  temperature  of  the  propellant  gases 

vQ  muzzle  velocity  of  the  projectile 

v  velocity  of  the  propellant  gases 

V  specific  volume  of  the  propellant  gases 


Ihm 

°R 


ft /sec 


p  pressure  in  the  gas  psia 

t  time  secs 


f0  time  at  start  of  gas  ejection 

B  force  on  breech 

a  speed  of  sound  in  gas 

tf«Cp/Cv  ratio  of  specific  heats  of  propellant  gas 

R  gas  constant 


secs 

lbf 

ft/sec 


lbf  ft 


1 


Symbol  Meaning 


Dimension 


g 


E 

KEP 

KEG 

KER 

KEB 

% 


g 


specific  heat  of  gas  at  constant 
pressure 

specific  heat  of  gas  at  constant 
volume 

Joule's  constant  «  777,5 


acceleration  due  to  gravity  «  32 , 17 
heat  of  explosion  of  the  propellant 


Btu 

tbsts 

Btu 


Item 

lbf  ft 


heat  of  explosion  of  the  igniter- 


energy  involved  in  a  process 

kinetic  energy  of  the  projectile  at  the 
start  of  gas  ejection  (tQ) 

kinetic  energy  of  the  gas  at  t0 

kinetic  energy  of  the  recoiling  parts 
at  tQ 

kinetic  energy  spent  in  engraving  the 
bore  and  in  forcing  the  projectile  thru* 
evaluated  at  t0 

thermal  energy  remaining  in  the  gas  at  t 


total  heat  released  by  charge  and  igniter 
during  combustion 

heat  released  by  charge 
heat  released  by  igniter 


lav 


dimensionless  velocity 


<6a  E _ 

T  pb 

V 

1"  vT 


init 

dimensionless  pressure  at  t0 
dimensionless  specific  volume  at  t 


Btu 


o 


Btu 

Btu 

^tem 

Btu 

Btu 

Btu 

Btu 

Btu 

Btu 

Btu 

BtU 

Btu 


o 
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Symbol 
&  _  T 


M 

P 

H 


J  = 


Meaning 

dimensionless  temperature  at  t0 ;  also 
used  as  an  angle  in  the  brake  design 

cumulative  mass  discharged 
cumulative  momentum  discharged 
cumulative  stagnation  enthalpy 


force  on  brake  in  axial  direction  of 
gun  tube 

vertical  force  on  brake  normal  to  axis 
of  gun  tube 

A  g  p 

- — -  — M11^  >  dimensional  constant  used 
b  init  “T  in  obtaining  a  dimension¬ 
less  time 

t  dimensionless  time 


Pg 


Dimension 

Ibm 

lbf  secs 
Btu 
Ibj 
Ibf 


—  ;  dimensionless  cumulative 
lb  init  momentum  discharged  from 
muzzle 


»  J®  dimensionless  cumulative  mass  discharged 
Mt 


sec 


-1 


dimensionless  cumulative  stagnation 
enthalpy  discharged 

axial  momentum  index,  i .  e . ,  the  ratio  of 
the  axial  brake  force  to  the  breech  force 

ratio  of  normal  brake  force  to  momentum 
rate  of  discharge 

-  1.46881 (1-  )  ratio  of  momentum 

propagated  in  axial  direction  thru  a 
control  surface  surrounding  brake  to 
total  momentum  discharged  from  the 
muzzle 


It 

b 


CO  : 


H 


B 

F 

f 


SYMBOLS  PERTAINING  TO  SHOCK  ESTABLISHMENT  AND  DECAY 


static  pressure  in  gas 

psia 

dynamic  pressure  in  gas 

psia 

specific  volume  of  gas 

ft3 

mr 

3 


Symbol  Meaning 


Dimension 


v 

v 

T 


*1  * 


V° 


!av 


/V 

r 


2 
M* 

r 

<*1 


velocity  of  gas 

velocity  of  center  of  mass  of  gas  mass 

static  temperature  of  gas 

Stagnation  temperature  at  center  of  shock 


specific  volume  at  center  of  shock  sphere 
^2av 

yo 

radius  of  shock  being  driven  by  mass  M 
discharged  from  weapon  (at  j”  =  i  ) 

total  mass  of  gas  within  sphere  of  radius 
r2 

radius  of  a  sphere  of  gas 

distance  from  center  of  muzzle  brake  side 
port  to  center  of  nearest  spherical  shock, 
at  maximum  shock  strength 

path  length  of  gas  flow  from  center  of 
muzzle  to  center  of  one  of  the  side  muzzle 
brake  ports 

horizontal  coordinate  axis  normal  to  gun 
tube  axis  with  center  at  center  of  front 
port  of  brake 

axis  normal  to  axis  of  gun  tube  and  X-axis, 
directed  up,  with  origin  at  origin  of 
X-axis 

coordinate  axis  along  axis  of  gun  tube 
directed  from  the  breech  with  origin  at 
center  of  front  port  of  muzzle  brake 

X-distance  to  center  of  a  side  port  of 
the  muzzle  brake 

Y-distance  to  center  of  a  side  port  of 
the  muzzle  brake 

Z-distance  to  center  of  a  side  port  of 
the  muzzle  brake 


"R 


'R 


31a 
ft  3/ lb 


m 


ft 

lb 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


m 
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Symbol  Meaning  Dimension 


XI 

X-distance  between  center  of  a  side 
port  of  m.  b.  and  center  of  associated, 
fully  developed  shock  sphere 

ft 

yi 

Y-distance  between  center  of  a  side  port 
of  m.  b.  and  center  of  associated,  fully 
developed  shock  sphere 

ft 

H 

Z-distance  between  center  of  a  side  port 
of  in.  b.  and  center  of  associated,  fully 
developed  shock  sphere 

ft 

X  =' 

x0+xl 

X-distance  to  center  of  shock  sphere 

ft 

y  = 

yo+yi 

Y-distance  to  center  of  shock  sphere 

ft 

z  =■ 

ZQ+Z1 

Z-distance  to  center  of  shock  sphere 

ft 

f 

horizontal  coordinate  normal  to  vertical 
plane  thru  gun  tube  axis  with  origin  at 
center  of  rear  trunnion 

ft 

vertical  coordinate  measured  positively 
upwards  with  origin  at  center  of  rear 
trunnion 

ft 

horizontal  coordinate,  in  plane  parallel 
to  ground,  orthogonal  to  fe  and  h  measured 
positively  toward  the  muzzle 

ft 

h 

height  above  £  -  »  plane  to  reference 
position  in  blast  field 

ft 

f 

distance  behind  center  of  rear  trunnion 
in  p  ->  V  plane  to  reference  position  in 
blast  field 

ft 

G 

distance  from  center  of  front  port  of 
muzzle  brake  to  center  of  rear  trunnion 

ft 

L 

distance  from  center  of  front  port  of 
muzzle  brake  to  reference  position  in 
blast  field 

ft 

N 

distance  from  center  of  nearest  shock 
sphere  at  max  static  overpressure  to 
reference  position  in  blast  field 

ft 

0 

distance  from  center  of  farthest  shock 
sphere  at  max  static  overpressure  to 
reference  position  in  blast  field 

ft 

D 

projection  in^-y  plane  of  distance  from 
center  of  front  port  of  muzzle  brake  to 
reference  position  in  blast  field 

ft 
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Symbol 


$ 


A 

ex' 


£ 


cos  A 


y 


cosA. 


cos  A. 


cosB 


f 


cosB, 


cosB, 


rad  or  deg 


Meaning  Dimension 

quadrant  elevation 


angle  between  gun  tube  axis  and  pro¬ 
jection  of  D  on  plane  of  tube  in  plane 
of  tube 

projection  of  in  Jt  -i>  plane 

angular  projection  in  X-Z  plane  that 
mean  flow  from  the  brake  makes  with 
Z-axis 

angular  projection  in  X-Y  plane  that 
mean  flow  makes  with  Y-axis 


direction  cosine  of  vj  vector  of  nearest 
shock  sphere  with  respect  to  t  -axis  in 
p  -  -  v  system  ' 

direction  cosine  of  v^  vector  of  nearest 
shock  sphere  with  respect  to  fj  -axis  in 
v  system 

direction  cosine  of  v^  vector  of  nearest 
shock  sphere  with  respect  to  P  -axis  in 
£  q  *  P  system 

direction  cosine  of  v^  vector  of  farthest 
shock  sphere  with  respect  to  £  -axis  in 
^  -  rj  -  p  system  > 

as  cosA  n  above  except  pertains  to 
farthest  shock  sphere 

as  cpsAy  above  except  pertains  to 
farthest  shock  sphere 


Pg  =  P“fto  static  overpressure  (subscripted) 


cp  *  Ps 
P 

K 


PF 

EFF 


dimensionless  static  overpressure 
(subscripted) 

effective  per  unit  Vj  contribution  to 
the  gas  velocity  behind  the  shock  front 

an  empirical  propagation  factor  defined 
in  equation  0,63 

per  unit  total  stagnation  enthalpy  in¬ 
volved  in  a  point  source  explosion  at 
the  muzzle;  efficiency  of  energy 
utilization  in  producing  shock 


psi 


psi 
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Symbol 

°T 

K  "  X 

DPP 

SPLSP 


DPLSP 


Meaning 

dimensional  constant  used  in  point 
source  explosion  theory 

non-dimensional  distance  from  center 
Of  point  source  explosion 

duration  of  the  positive  static 
pressure  phase  of  the  shock  wave 

impulse  per  unit  area  from  the 
positive  static  pressure  phase  of 
the  shock  wave  due  to  a  point 
source  explosion 

impulse  per  unit  area  from  the  posi¬ 
tive  dynamic  pressure  phase  of  the 
shock  wave  due  to  a  point  source 
explosion 


Dimension 

ft 

millisecs 

psi  millisecs 

psi  millisecs 
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X  :>*  N 


SUBSCRIPTS 


o  refers  to  conditions  at  the  muzzle 

oo  refers  to  ambient  (free  field)  conditions 

g  pertains  to  the  combustion  products 

b  refers  to  value  of  variable  at  the  breech 

eff  refers  to  an  effective  value 

av  refers  to  an  average  (or  weighted  average) 

value  over  a  volume 

init  refers  to  initial  conditions  for  gas  ejection 

T  refers  to  a  cumulative  or  integrated  value 

x  refers  to  conditions  upstream  of  a  normal  shock 

y  refers  to  conditions  downstream  of  (or  behind) 

a  normal  shock 

pertains  to  components  along  X-,  Y-,  and  Z-, 
coordinate  axes 


0  refers  to  location  of  muzzle  brake  ports  in 

X-Y-Z  system 

1  pertains  to  each  of  the  two  spherical  shocks 
created  by  the  muzzle  brake,  at  max  peak 
overpressure 

2  pertains  to  conditions  in  a  shell  of  compressed 
ambient  air  in  the  single  shock  sphere  model  as 
observed  in  a  frame  stationary  with  respect  to 
the  center  of  mass  of  this  system 

3  refers  to  conditions  at  reference  position  in 
blast  field 
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SUPERSCRIPTS 


refers  to  a  stagnation  value 

distinguishes  angular  projections,  as  noted 

refers  to  conditions  predicted  by  the  point 
source  explosion  model 

pertains  to  motion  of  center  of  mass  of  gas 
discharged 

refers  to  normal  shock  conditions  as  viewed 
by  an  observer  with  respect  to  whom  the  shock 
is  moving  in  still  air 


INTRODUCTION 


It  has  appeared  to  be  desirable  to  have  an  analytical 
means  for  predicting  the  severity  of  the  artillery  blast 
within  the  crew  area.  The  practice  of  designing  muzzle 
brakes  from  only  the  point-of-view  of  reducing  rod  pull 
and  without  regard  to  the  effect  of  the  brake  design  upon 
the  blast  field  has  been  perpetuated  because  a  good 
analytical  tool  for  predicting  the  performance  of  a  brake - 
weapon  system  in  terms  of  operator  comfort  has  been  lack¬ 
ing.  The  purpose  of  this  report  is  to  fill  that  need 
insofar  as  engineering  results  are  concerned.  There  is 
much  more  that  one  could  say  that  would  give  a  better 
understanding  of  the  physical  mechanisms  involved  in  the 
creation  of  atmospheric  shock  from  artillery. 

It  is  recognized  that  the  mathematical  model  presented 
here  has  oversimplified  the  phenomenon.  At  the  outset  of 
this  study  there  was  the  hope  that  the  assumptions  involved 
in  the  model  were  not  so  unrealistic  that  quantitative 
agreement  with  experiment  was  precluded.  As  a  working 
goal;  it  was  decided  that  an  agreement  in  the  peak  static 
overpressure  of  better  than  20%  of  max  value  for  all 
points  within  the  crew  area  was  necessary.  To  achieve 
this  goal,  a  minimum  number  of  empirical  constants  were 
to  be  introduced  and  ad  hoc  assumptions  were  to  be  elim¬ 
inated,  if  possible.  . 

Since  there  seems  to  be  a  positive  correlation 
between  the  discomfort  that  a  human  operator  experiences 
and  the  peak  static  overpressure  in  the  incident  shock 
wave,  and  since  the  latter  parameter  has  been  accepted 
for  some  time  as  a  critical  parameter  in  the  man-weapon 
system,  this  theory  was  constructed  to  predict  the  peak 
static  overpressure  for  any  position  within  the  crew  area. 
Additionally,  the  peak  dynamic  pressure  component  is 
given. 

For  reference  and  comparison  purposes  the  following 
parameters  predicted  by  a  competitive  theory  are  given: 
the  peak  static  and  dynamic  pressures,  the  total  positive 
impulse,  and  the  duration  of  the  positive  pressure  phase 
in  the  shock  wave.  These  computations  were  made  by 
assuming  a  point  source  explosion  at  the  muzzle  with  a 
total  energy  release  equal  to  the  total  stagnation 
enthalpy  discharged  times  an  efficiency  factor .  The 
latter  values  indicate  what  is  expected  for  a  spher¬ 
ically  symmetric  burst  and,  by  comparison  with  the 
results  of  the  present  theory,  indicate  the  directional 
effects  of  the  brake-weapon  system. 
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In  predicting  the  peak  static  overpressure,  it  was 
desirable  to  take  into  account  specific  geometric  details 
of  the  muzzle  brake  such  as  the  size  and  shape  of  the 
brake  and  to  take  into  account  such  physical  factors  as 
the  momentum  indices  in  the  axial  and  transverse  vertical 
directions  (for  vertically  asymmetric  brakes).  This  was 
done  in  sufficient  generality  to  include  the  majority  of 
the  presently  used  brake  designs.  Also,  the  interior 
ballistic  factors  were  considered  as  well  as  the  relevant 
weapon  dimensions. 

Computations  were  made  from  the  theory  using  a  digital 
computer.  The  computed  results  have  been  compared  (insofar 
as  is  possible)  with  recent  experimental  results  obtained 
by  the  Artillery  weapons  Branch,  Artillery  Division, 
Development  and  Proof  Services  for  the  Human  Engineering 
Laboratories,  Aberdeen  Proving  Ground,  Md.  (Reference  12.) 

Since  hand  calculation  from  the  theory  is  extremely 
laborious,  the  author  has  included  the  logic  diagram  and 
the  FORTRAN  program  for  making  these  computations  with  a 
digital  computer . 


THE  MATHEMATICAL  MODEL 


It  is  recognized  that  gaseous  discharge  from  the 
muzzle  of  an  artillery  piece  does  not  admit  the  same 
simplifications  in  the  mathematical  model  as,  say,  the 
explosion  at  a  point,  in  which  case  an  instantaneous 
release  of  energy  is  validly  assumed.  The  discharge 
from  a  gun  tube  is  neither  a  pure  energy-  nor  a  pure 
mass-source  as  far  as  the  shock  field  near  the  muzzle 
is  concerned. 

In  the  case  of  gas  discharged  from  a  gun  tube 
without  a  muzzle  brake,  the  events  occurring  prior  to 
the  establishment  of  maximum  shock  strength  have  been 
described  in  some  detail  qualitatively  in  reference  7. 
The  essential  feature  of  that  description  which  I  wish 
to  use  is  the  idea  of  a  "shock  bottle. "  As  the 
combustion  products  are  discharged,  they  push  back  the 
ambient  air  replacing  a  volume  formerly  occupied  by 
air.  This  process  may  be  thought  of  as  a  three  dimen¬ 
sional  fluid  piston  expanding  in  air.  During  this 
discharge,  which  occurs  transonically  in  the  combustion 
products  but  supersonically  relative  to  the  air,  a 
shock  wave  is  established  in  the  air.  This  shock 
grows  in  strength  as  more  gas  is  fed  into  the  bottle¬ 
shaped  volume  of  combustion  products  at  the  muzzle . 
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F  I  G  U  R  E  1.  A 

S  I  N  G  L  B  S  P  H  E  R  E  S  M  O  C  K  M  O  D  E  L 


A DVANC I  N G  S  NOG  K 


1 .  B 


Since  the  mass-*  momentum-  and  enthalphy-rates  of 
discharge  are  rapidly  decreasing,  a  point  is  reached 
where  the  bottle  of  driving  gases  cannot  affect  an  in¬ 
crease  in  shock  strength,  i . e. ,  cannot  accelerate  the 
shell  of  compressed  air  advancing  before  it.  At  this 
moment  in  time  the  shock  is  said  to  be  fully  developed. 
Subsequent  propagation  of  the  shock  wave  in  the  air 
away  from  the  driving  gases  only  serves  to  attenuate  the 
peak  overpressure,  baring  coalescence  of  two  shock  waves 
upon  shock  reflection. 

At  maximum  shock  strength  the  dimensions  of  the 
shock  bottle  are  such  that  the  max  streamwise  dimension 
from  the  muzzle  is  approximately  equal  to  the  max  trans¬ 
verse  dimension  (reference  7) . 

For  weapons  such  as  the  standard  105mm  howitzer, 
this  has  been  found  to  occur  at  a  non-dimensional  time 
of  -  1 .  Certainly,  one  would  suspect  that  this  time 
would  vary  from  weapon  to  weapon  and  from  one  zone 
charge  to  another,  for  a  given  weapon .  However,  slight 
variations  in  the  predicted  time  of  max  shock  strength 
should  not  alter  appreciably  the  predicted  max  over¬ 
pressure,  since  the  latter  value  changes  slowly  in  time 
hear  its  maximum,  where  time  derivatives  are  zero* 

Thus,  in  the  mathematical  model,  max  peak  over¬ 
pressure  occurs  at  7“**  1.  Further  simplification  in 
the  model  was  obtained  by  assuming  the  shock  bottle  is 
spherical.  The  geometry  of  the  shock  model  is  illustrated 
in  figure  la. 

Variables  within  the  driving  gas  at  the  interface 
with  the  shell  of  compressed  air  are  subscripted  with 
the  number  2.  Conditions  within  the  thin  Shell  of 
compressed  air  are  assumed  to  be  uniform  and  are  sub” 
scripted  with  the  letter  y .  Conditions  upstream  of  the 
moving  shock  ambient  conditions,  in  this  case  — 
are  subscripted  with  the  letter  x  or  the  symbol  . 

Since  the  combustion  gases  are  also  transporting 
momentum  along  the  tube  axis,  the  shock  sphere  must  be 
moving  in  an  axial  direction  with  a  mean  velocity  -‘Vg  • 

In  a  frame  of  reference  moving  with  the  center  of  mass, 
the  state  variables  are  assumed  to  have  the  distribu¬ 
tions  shown  in  figure  lb. 

Note  Py  ~  p2,  i . e . ,  that  at  max  shock  strength  the 
pressure  gradient  at  the  driving  interface  has  vanished. 
Also  note  that  turbulent  mixing  has  not  been  such  as  to 
raise  the  temperature  of  the  compressed  air  shell  to 
the  value  T2.  This  assumption  is  based  upon  the  fact 
that  maximum  shock  strength  is  developed  in  a  very  brief 
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interval  and  that  heat  transfer  mechanisms  in  a  dilute  gas 
require  a  longer  active  period  to  be  effective.  A  third 
significant  aspect  of  this  model  is  the  requirement  that 
the  mach  number  of  the  driving  gas  at  the  interface  equal 
unity.  If  the  mach  number  were  lower,  expansion  would 
quickly  raise  its  value ;  however,  if  the  mach  number  here 
were  higher  than  unity,  shock  reflection  at  the  interface 
would  bring  its  value  down. 

The  value  of  the  stagnation  pressure  and  of  the 
stagnation  temperature  of  the  combustion  gases  within  the 
driving  sphere  and,  in  fact,  the  size  of  the  sphere  it¬ 
self  is  a  function  of  the  cumulative  -mass,  -momentum, 
and  -stagnation  enthalpy  discharged  up  to  S’  =  1.  Further 
the  Rankine-Hugoniot  conditions  (Appendix  I)  require  that 
definite  relationships  exist  between  variables  on  either 
side  of  the  shock  —  assumed  here  to  be  a  one-dimensional 
normal  shock  —  e . g. ,  between  py  and  px.  Conditions  in 

front  of  the  shock  are  known.  Therefore,  one  can  write  a 
set  of  equations  which  determine  the  peak  overpressure, 
Py-Px  *  Psi  at  maximum  shock  strength.  The  derivation  of 
these  equations  is  found  in  Appendix  II .  It  is  also  shown 
in  this  appendix  that  the  presence  of  a  muzzle  brake,  which 
alters  the  distribution  of  mass  from  the  muzzle,  does  not 
change  the  maximum  peak  overpressure  in  the  stationary 
shock  frame  of  reference.  It  is  assumed  that  the  muzzle 
brake  causes  two  initially  distinct  shock  spheres  to  be 
developed,  each  of  which  propagates  with  an  independent 
decay  rate  until  dpalescence.  Constructive  interference 
within  the  crew  area  is  assumed .  (See  Appendix  HI.) 

The  effect  of  shock  reflection  from  the  ground  is  not  con¬ 
sidered,  altho  at  low  quadrant  elevations  this  phenomenon 
undoubtedly  plays  a  role  in  determining  peak  overpressure 
in  the  crew  area. 

The  directional  effects  of  the  shock  wave  are  ob¬ 
tained  by  considering  the  maximum  peak  overpressure  in 
a  frame  of  reference  that  is  fixed  with  respect  to  the 
ground .  Reference  3  suggests  that,  for  the  simplified 
model  considered  there,  strong  directional  effects  are 
absent  in  the  blast  field  of  a  moving  fluid  source  when¬ 
ever  the  mach  number  of  the  source  in  the  medium  is  less 
than  unity.  However,  even  for  a  subsonic  source,  there 
are  directional  effects.  These  result  in  a  strengthening 
of  the  shock  in  the  direction  of  motion  of  the  source  and 
a  weakening  of  the  shock  in  the  opposite  direction.  To 
account  for  the  directional  effects  due  to  the  motion  of 
the  shock  sphere  relative  to  the  ground,  the  center-of- 
mass  velocity  component  in  the  direction  of  the  reference 
position  was  multiplied  by  an  empirical  propagation 
factor  and  added  to  the  max  radial  velocity  of  the  shock 
wave  in  the  center-of-mass  system  to  obtain  a  maximum 


effective  shock  wave  velocity  in  a  ground  fixed  frame , 

By  the  Rankine-Hugoniot  relations,  the  peak  overpressure 
is  dependent  upon  the  shock  velocity.  Thus,  the  maximum 
peak  overpressure  in  the  ground  frame  of  reference  can 
be  obtained. 

Having  obtained  the  max  peak  overpressure,  one  must 
consider  the  rate  of  decay  of  a  shock  wave  with  distance 
to  be  able  to  predict  the  peak  overpressure  at  the 
reference  position.  Reference  3  indicates  that  over¬ 
pressure  varies  approximately  as  r”1  as  the  overpressure 
approaches  zero.  Reference  1,  considering  point  source 
explosions  in  air,  indicates  a  much  more  rapid  decay  for 
very  strong  shocks  with  an  approximately  r“2/2  dependence 
for  shocks  having  a  strength  within  the  range  commonly 
encountered  in  muzzle  blast  (7  -  .1  atm  overpressure) , 
the  dependence  gradually  approaching  r”*  for  shocks  of 
infinitesimal  strength . 

It  seems  appropriate  to  choose  an  r“2/2  overpressure 
decay  for  the  present  mathematical  model .  Then,  knowing 
the  max  peak  overpressure  and  the  shock  radius  at  this 
time  and  knowing  the  distance  between  the  center  of  the 
shock  and  the  reference  position,  one  can  compute  the 
peak  overpressure  at  the  reference  position.  Details  of 
the  derivation  of  the  equations  concerning  the  over¬ 
pressure  at  the  reference  position  are  found  in  Appendix 
III. 

For  the  derivation  of  the  equations  concerning  the 
gun  internal  ballistics,  the  reader  is  referred  to  the 
appendicies  of  reference  10. 

The  equations  referring  to  conditions  in  the  blast 
field  due  to  a  point  source  explosion  were  obtained  from 
references  1  and  2.  It  was  found  that  quadratic  best 
fits  on  log-log  paper  to  the  data  given  in  these  refer¬ 
ences  provided  convenient  computational  formulae,  of 
sufficient  accuracy,  for  obtaining  the  desired  parameters . 
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EQUATIONS 


The  following  equations,  used  in  the  computer  program, 
(listed  in  order  used)  are  derived  in  Appendices  I,  II, 
and  III  of  this  report  and  in  Appendices  I,  II,  and  III  of 
reference  10. 
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THE  DIGITAL  COMPUTER  PROGRAM 


The  computer  program*  designed  to  evaluate  the  equa¬ 
tions,  permits  one  to  choose  the  parameter  space  for  which 
the  computations  of  pressure,  etc .  are  made.  One  must 
enter  the  relevant  gun  constants  and  internal  ballistic 
parameters  (listed  later  in  this  section)  along  with  the 
desired  initial  values  and  increments  of  the  following 
variables:  to  ,  V  ,  QE,  h,  f ,  tS-  »  Provision  has  been 
made  for  specifying  the  number  of  values  of  each  of  the 
above  variables .  The  number  of  values  (levels)  for  each 
variable  is  given  a  code  word  in  the  FORTRAN  program. 

The  values  are  called  in  order  into  computer  memory  under 
the  following  code  words: 


Variable 

Code  Word 
for  Counter 

Code  Word  for 
Number  of  Levels 

OS) 

KATE 

KENDAL 

V 

KAY 

KEN 

QE 

NIEL 

NORMA 

h 

MAC 

MAGE 

f 

LENNY 

LAURA 

# 

IRENE 

IVAN 

Thus, 

if  one  required  answers 

for  two  levels  of  to  , 

would  enter  2  for  the  code  word  KENDAL ,  Operationally,  the 
computer  would  establish  a  counter  labeled  KATE  which  would 
take  on  integral  values  from  1  thru  KENDAL,  or  2,  in  this 
case . 

The  computer  performs  a  complete  factorial  among  the 
variables.  Variables  are  listed  in  the  above  table  in 
order  of  progressive  nesting,  *$-  being  most  deeply  nested. 

Using  the  notation  Ax  as  the  finite  increment  of  the 
variable  x,  one  must  prepare  data  cards,  in  the  following 
order,  having  the  values  of :  C  ,  h,  f,  -fb  ,  QE,  V  ,  Ah, 

Af,  A ,  AQE,  A  P  ,  a,  G,  xp,  Yos  Zp,  co  ,  Poo  *  PWR,  A  to  , 

PF,  EFF,  where  PWR  represents  the  abs  value  of  the  power 
of  distance  in  the  peak  overpressure-distance  relationship 
—  taken  here  as  3/2.  The  tables  of  X^  (T°) ,  Vy/Vx  (T°) ,  Pa(l^) 
follow  the  above  data.  Finally,  the  gun  constants  and 
internal  ballistic  parameters  are  entered  as  follows; 

*2*  R2-  ec"  eig^  A?  Acnct-"  ^T*  Mc^  MigJ  MpF  kIR'  vo  " 

A  logic  diagram  of  the  program  and  the  source  program, 
written  in  full  FORTRAN,  follow. 
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NUMBER  STATEMENT 

PEAK  OVERPRESSURE  SCHLENKER-OLSON  RIA  R+D  OCTOBER  1962 
READ  INPUT  TAPE  1, 819,  KENDAL, KEN, NORMA,  MAGE, LAURA, IVAN 
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FH2=EMSTR*VlBAR**2/(2. *32. 17) +C0EF2*T0#144. *PINE*C0EF1*CYE1*(1. +VYVX) 


IF  (FTO) 18, 19, 19 


o 

Eh 

Ixj 

I 

s— / 
<1 
e 


00 
I— ( 


rH 

o 

Eh 


o 

H 

Ph 

& 

Eh 

Ph 


© 

CM 

,<»> 

O 

© 


o 

Q 


o 

Eh 

pH 

03 

< 


m 


CO 

CM 

*\ 

CM 

CM 

«\ 

rH 

CM 


O 

Eh 


0> 


© 

CM 


o 

Eh 

l 


rH 

CM 


O 

Eh 

Ph 


5  8  5  5  8  S  £ 


O 

Eh 

Ph 

CQ 

< 

A 

S3 


Eh 

Q 


O 

Eh 


o 

Eh 


CO 

IL 


o 

o 

Eh 

I 


O 

O 

CO 

Eh 

CM 

w 

rH 

\ 

** 

/■N 

© 

o 

co 

o 

CM 

Eh 

■rH 

pH 

1 

iO 

o 

CM 

Eh 

1*4 

o 

Eh 

1*4 

' — ' 

* 

\ 

o 

Q 

o 

CO 

5 

Eh 

o 

Eh 

cO 

Eh 

o 

1*4 

Eh 

O 

P*4 

Eh 

1*4 

O 

O 

V 

H 

o 

1 

EH 

Eh 

CM 

EH 

Eh 

6 

o 

II 

II 

O 

n 

Q 

Q 

1*4 

oq 

G 

o 

Eh 

O 

Eh 

Q 

H4 

M 

Eh 

Eh 

Ph 

CO 

m 

CM 

CM 

CM 

CM 

Ch 

rH 


32 


ABF=ABFT0 


i 

« 

Er»- 


Eh 


A 

00 

/A 

\ 

rH 

• 

1 

rH 

>Q 

w 

W 

* 

W 

* 

Eh 

As 

l 

An 

X 

n 

|jj 

u 

Eh 

M 

© 

Cr 

Or 

Eh 

* 

* 

•> 

CD 

/Hn 

• 

CM 

CM 

r\ 

rt> 

CQ 

rH 

rH 

rH 

rfi 

Or 

O 

rH 

n 

HJ* 

•rH: 

>0 

© 

(0 

•«* 

N-r 

\ 

© 

CM 

03 

CQ 

in 

PO 

< 

rH 

© 

ft 

'• 

•* 

M 

•> 

rH 

I 

fir 

!> 

c*> 

•V 

/H 

fi"“ 

/TN 

CO 

03 

*"5 

fi- 

Pa 

0? 

© 

* 

ft 

<o 

rH 

rH 

CQ 

CQ 

53 

03 

/T"\ 

ft 

< 

Eh 

rH 

o 

© 

/•> 

rH 

o 

•v 

EH 

rr\ 

rH 

Hh 

00 

pa 

Eh 

rH. 

[*r 

Eh 

% 

w 

/■■n 

CO 

* 

P 

CO 

m 

Eh 

P 

CM 

fc< 

tH 

CO 

Ph 

ft 

03 

< 

O 

P 

O 

W 

CO 

I 

H 

1 

•M* 

pH 

EH 

rH 

1 

O 

EH 

w 

CQ 

a 

CM 

►<? 

1 

M 

P 

• 

pa 

p 

s-* 

pL4 

Eh 

fa 

Q 

CQ 

i 

rH 

& 

P 

rH 

S-rr 

o 

CQ 

o 

h 

m 

+ 

w 

Q 

rH 

O 

!5 

03 

\ 

n 

© 

< 

Eh 

n 

<2 

O 

o 

J 

o 

Eh 

o 

H 

t-H 

Or 

(H 

< 

pa 

H 

03 

1 

6 

g 

Eh 

Eh 

EH 

Eh 

B 

rH 

s-< 

g 

Or 

X 

Eh 

1 

6 

o 

pH 

I 

p 

« 

« 

04 

A 

ft 

tH 

O 

< 

pL< 

o 

EH 

CQ 

O 

Q 

w 

q 

O 

© 

Pm 

o 

03 

□ 

« 

O 

HH 

W 

Eh 

fa 

Eh 

O 

p 

a 

Eh 

a 

M 

Q 

Or 

w 

< 

© 

rH 

© 

<0 

© 

CM 

03 1 

CM 

rH 

rH 

tH 

O 

© 

o 

CM 

rH1 

rH 

CM 

O 


33 


800  FORMAT  (3HPS2, 4X, F7 . 2, 5H  PS IG, 1 IX, 2HTO, F8 . 2,  2H  R///) 
R2TIL-3 . *C0EF1 *C YE 1/ ( 4 . *3. 1415927) * (1 . +VYVX) 

Rl»  (R2TIL/2..  )  **  ( 1 .  /3 .  ) 


R2TIL=R2TIL**(1 ./8. ) 
V2-A*SQRTF (6 . *PS2/ (7 . +PINE )+l. 
0MEGA=A2 


% 


s. 

• 

00 

a 

• 

ft 

CO 

* 

ft 

H 

«V 

a 

< 

X 

H 

P 

o 

m 

< 

rH 

B 

o 

ft 

s. 

s 

rH 

«G 

ft 

Eh 

CM 

Q 

\ 

Eh 

CQ 

« 

* 

hg 

ft 

9\ 

ft 

* 

o 

© 

N 

Eh 

P 

ft 

rH 

rH 

rH 

cm 

ft 

1 

< 

i 

i 

P 

1 

a 

a 

n 

00 

a 

CO 

ft 

ft 

ffl 

a 

a 

> 

o 

CM 

ft 

1 

rH 

a 

a 

a 

1 

■ft 

1  ^ 

a 

in 

CQ  H 

iH 

ft 

ft 

a 

CM 

ft 

ft 

CM 

< 

m 

• 

1 

a 

a 

* 

»k 

© 

• 

* 

a 

1  w 

* 

•  a 

a 

rH 

» 

X 

a 

rH 

rH 

* 

■ts 

i 

X! 

- — s 

P!“ 

i  a 

a 

ft 

a 

?> 

\ 

\ 

ft 

ft 

in 

ft 

ft 

ft 

1 

o 

CM 

a 

BS 

ft 

s 

!  Eh 

p 

a 

a 

* 

0 

. 

rH 

ft 

rH' 

P 

hS 

<13 

ft 

ft' 

ft 

ft 

rH  rH 

rH 

N 

Ss 

a 

a 

P 

1 

,H 

I  Eh 

\ 

ft 

K 

a  \ 

\ 

* 

Eh 

\ 

/S 

ft 

ft 

CM 

II 

to 

\ 

< 

00  < 

< 

rH 

P 

\ 

rH 

rH 

> 

> 

* 

W  ft 

N 

\ 

rH 

CO  P 

o 

N 

ft 

\ 

o 

rH 

X 

* 

* 

* 

E- 

!  e 

\ 

fa 

ji 

ft  a 

ft 

1 

e 

v< 

00 

rH 

SH- 

< 

«3 

a 

< 

5 

v 

Jz; 

rH 

ft 

<  ft 

a 

ft 

p 

<C!\ 

rH 

ft 

§ 

o 

««3 

X 

i  O' 

M 

<$ 

< 

r— 

I  <5 

o 

P 

o 

rH 

rH 

Eh 

a 

Eh 

CO 

11 

» 

ii 

* 

* 

* 

Eh 

X 

rH 

a 

>4 

psj 

rH 

Q 

w 

< 

< 

ft 

< 

!  rH 

ft 

■rH 

ft 

v 

O’ 

< 

o 

ft 

rH 

1  Eh 

a 

w 

<3 

e 

!  ° 

© 

P 

Eh 

A 

a 

ii 

II 

|| 

M 

a 

H 

ft 

s 

i  II 

l| 

|| 

M 

a 

O 

n 

ft 

ft 

ft 

S 

1  tf 

o 

rH 

ft 

2 

>— 

rH 

ft 

ft 

a 

o 

Ph 

Eh 

ft 

Csl 

fen 

X 

a 

1  ^ 

ft 

Q 

a 

p 

< 

I  CSI 

ft 

X 

Se 

ft 

M 

CO 

X 

> 

ft 

o 

p 

■rH- 

CM 

CM 

p 

o 

O 

a 

ft 

a 

a 

34 


© 

© 

X 

O 

00 

•» 

CO 

© 

00 


pH 

M 


N 

> 


X 

pi 


04 

04 

0) 

04 

04 

Q, 

5 

H 

CO 


X 

> 

B 

8 

M 

« 

rH 

X: 

X 


X 

© 

00 

•l- 

© 

o 

© 

•* 

© 

© 

Q0 

i^S 

rH 

N 

X 


X 

+ 

X 

<N 

© 

in 

rl 

X 


CO 

l 


oo 

© 

GO 

O 


04 

s. 

X 

04 

© 

© 

rH 

X» 

rH 

• 

Q0 


m  ^  8  i 


© 

© 

rH 

© 

/-N 

X 

C3 

© 

rH 

© 

© 

X 

•i 

X 

© 

X 

>s 

© 

X 

04 

04 

© 

rH 

rH 

\ 

s. 

• 

© 

x 

© 

rH 

X 

X 

04 

> 

Q| 

x 

04 

© 

• 

© 

X 

© 

* 

X 

© 

VwT 

© 

X 

Mi 

s 

oo 

pH 

rH 

X 

* 

O0 

X 

/■> 

z 

rH 

X 

a 

£ 

© 

EH 

rH 

X 

rH 

i — ( 

j 

X 

oo 

X 

X 

EH 

© 

• 

X 

a 

> 

X 

co 

i 

a 

o 

H 

1 

o 

« 

J 

EH 

X 

Eh 

Eh 

X 

X 

W 

x. 

X 

a 

X 

O 

3 

Ph 

X 

Q 

a 

A 

a 

© 

3 

w 

X 

o 

a 

X 

a 

CO 

o 

© 


© 

© 

X 

© 

© 

© 

© 

© 

© 

© 

© 

o 

© 

o 

© 

© 

© 

O  rH 
rH  rH 
©  00 


35 


WRITE  OUTPUT  TAPE  7,  812, ALHAT, BHAT, D1 
812  FORMAT  (60X,  9HALPHA-H AT ,  6X,  8HBETA-HAT,  7X,  2HD1/60X,  F7 2,  4H  PEG,  4X 
F6. 2j  4H  DEG,  5X,F6.2,3fl  FT////) 

QE-B1 

DO  6  NIEL“1, NORMA 


WRITE  OUTPUT  TAPE  7,813,QE 


* 

* 

oo 

p* 

in 

01 

eq 

t* 

© 

s 

< 

Eh 

1 

Eh 


© 

W 

i* 

§ 

M 


00 

e* 

in 

01 

CQ 


CM 

* 

★ 

do 

IS 

•n 

01 


vrs 

■fjs. 

fa 

lift 

© 

>pn 

w 

a 

UJ 

s 

pH 

o 

Q 

a 

© 

© 

w 

<y 

>S 

fa 

a 

/r< 

a 

w 

a 

© 

EH 

Eh 

HJ* 

fa 

Qf 

© 

© 

S3 

fa 

fa 

© 

w 

*S 

CM 

o 

fa 

• 

•> 

m 

a 

Eh 

• 

o 

© 

P* 

M 

CO 

© 

U 

© 

* 

o 

© 

• 

•» 

•» 

a 

fa 

o 

* 

(0 

S3 

CM 

X 

© 

?> 

© 

* 

>s 

fa 
-  •* 

•» 

• 

•> 

• 

•» 

X 

is 

N 

a 

© 

© 

© 

© 

in 

© 

* 

CM 

E«i 

rH 

rH 

fa 

pH 

© 

* 

* 

© 

© 

00 

* 

© 

•\ 

© 

\ 

CM 

/Hn 

* 

o 

01 

•» 

X 

* 

X 

<3 

• 

© 

a 

o 

\ 

p- 

in 

© 

fa 

.Eh- 

fa 

fa 

Eh 

\ 

p» 

< 

© 

a 

© 

© 

© 

a 

w 

m 

fa 

\ 

X 

X 

\ 

a 

a 

a 

\ 

© 

o 

eh 

O' 

X 

P 

fa 

fa 

<! 

fa 

h 

<■ 

CM 

© 

a 

X 

W 

S3 

X 

< 

< 

a 

> 

S3 

a 

Eh 

• 

1 

M 

CM 

| 

Eh 

rH 

J 

Eh 

pH 

M 

Eh 

© 

S 

fa 

© 

• 

© 

H 

•> 

pi 

•> 

ft 

a 

rH 

<s 

X 

a 

X 

rH 

Eh 

X 

.pH1 

Eh 

X 

Eh 

\ 

Si> 

fa 

© 

p 

X 

J 

P 

» 

P 

© 

>< 

< 

a 

Ol 

rH 

& 

t> 

fa 

& 

© 

fa 

«< 

fa 

Eh 

Eh 

< 

§2 

Ss 

Q 

Eh 

\ 

a 

H 

\ 

Eh 

\ 

© 

a 

a 

V 

P 

V 

X 

P 

s>/ 

a 

P 

w 

© 

a 

O’ 

«*! 

o 

a 

O 

a 

a 

Q 

O 

Eh 

© 

J 

Eh 

< 

a 

Eh 

J 

H 

I 

w 

£ 

R 

w 

a 

a 

Eh 

< 

a 

< 

< 

a 

§ 

a 

fa 

2S 

a 

<N 

•n 

a 

w 

© 

Eh 

a 

X 

Eh 

© 

Eh 

a 

a 

Eh 

© 

Eh 

a 

X 

ffl 

M 

es 

© 

*H 

a 

M 

© 

O 

a 

H 

o 

1 

o 

X 

o 

jf 

2 

X 

o 

a 

E-i 

2 

a 

o 

O 

u 

M 

s 

fa 

S3 

Q 

a 

fa 

fa 

Q 

a 

fa 

Q 

a 

fa 

O 

+ 

© 

Eh 

CO 

■o* 

© 

rH 

(— 1 

rH 

© 

© 

© 

© 

rH 

00 


36 


EH  (G*COSF (QE/57 . 29578) +F) /COSTH 


817  FORMAT  (45X,  9HTHET A-HAT ,  10X,  1HD/45X,  F6 . 2,  4H  DEG,  5X,  F6 . 2,  3H 


fa 


* 

eg 

+ 

fa 

EH 

0 

O 

0 

* 

(A 

PH 

VP 

+ 

CM 

# 

* 

pa 

fa 

■| 

A 

00 

IS 

SB 

<35 

CM 

IS 

CM 

A 

SB 

* 

CM 

\ 

* 

* 

m 

As 

* 

at 

A 

PA 

VP 

00 

fa 

fa 

P- 

i 

03 

SB 

A 

O 

0 

oo 

0 

CM 

s 

SB 

* 

• 

!s 

p* 

0 

4- 

SB 

CM 

A 

\ 

A 

• 

00 

fa 

A 

a 

P» 

s> 

at 

a 

o 

SB 

iB 

VP 

o 

0 

\ 

<35 

fa 

0 

4- 

fa 

CM 

X 

4- 

CM 

o> 

• 

M 

CM 

* 

■VP 

O 

0 

* 

* 

Cm 

0 

* 

•* 

PA 

a 

\ 

>< 

PA 

X 

PH 

fa 

1 

X 

1 

CO 

at 

PA 

4- 

fa 

* 

vp 

oo 

fa 

& 

Eh 

0 

fa 

55 

■w? 

a 

in 

£ 

PH 

4- 

< 

ph 

0 

w 

0 

8 

* 

a 

CO 

* 

CM 

• 

0 

* 

* 

8 

8 

J 

PA 

tv 

8 

vp 

1 — / 

eg 

SB 

v 

VP 

fa 

\ 

4* 

VP 

fa 

pH 

pH 

pH 

Ei 

p3 

0 

fa 

fa 

EH 

X 

>H 

(S3 

fa 

fa 

vp 

at 

Eh 

fa 

+ 

•¥ 

4> 

O’ 

13 

VP 

V 

fa 

at 

o 

© 

O 

03 

fa 

6 

fa 

cy 

0 

X 

is 

N 

is 

>» 

a 

0 

0 

1 

1 

H 

r 

-4 

< 

o 

O 

i 

a 

X 

t* 

S3 

w 

o 

0 

6 

fa 

p-i 

fa¬ 

O 

a 

's 

\ 

fa 

A 

o 

pA 

A 

\ 

A 

P”V 

0 

0 

fa 

IS 

a 

PH 

0 

■i 

0 

A 

0 

PA 

0 

0 

CM 

0 

CM 

A 

tA 

o 

l> 

b 

0 

0 

0 

PH 

0 

P? 

P“ 

P» 

05 

0, 

0 

0 

0 

0 

CM 

\ 

CM 

N. 

0 

<35 

a 

fa 

b 

fa 

CM 

0 

tH 

c? 

PH 

cy 

0 

VP 

0 

V 

IS 

PH 

\ 

fa 

\ 

fa 

0 

0 

fa 

a 

fa 

a 

\ 

\ 

cy 

M 

cy 

M 

fa 

fa 

VP 

0 

VP 

8 

cy 

at 

fa 

* 

fa 

* 

VP 

VP 

0 

pH 

0 

!h 

fa 

fa 

o 

1 

o 

! 

a 

a 

0 

PA 

0 

A 

PH 

PH 

* 

Q0 

# 

0 

8 

0 

i* 

tv 

>H 

PH 

* 

# 

4- 

0 

4* 

0 

pH 

pH 

A 

0 

A 

0 

eg 

is 

CD 

CM 

00 

CM 

> 

> 

tH 

0 

p- 

• 

4- 

3 

0 

0 

P* 

PA 

PA 

<35 

0 

0 

0 

0 

0 

CM 

Ss 

CM 

s 

IS 

IS 

0 

fa 

0 

w 

0 

0 

tv 

at 

P* 

cy 

0 

(35 

0 

VP 

0 

SP 

CM 

0 

\ 

fa 

\ 

fa 

c 

o 

a 

fa 

0 

fa 

0 

N 

tv 

h 

at 

o 

o 

at 

o 

0 

0 

\ 

VP 

0 

\ 

VP 

0 

\ 

\ 

PA 

fa 

* 

fa 

■* 

fa 

fa 

X 

a 

fa 

X 

a  i 

P3 

cy 

cy 

PH 

4? 

4h 

M 

ir- 

VP 

V 

fa 

0 

fa 

fa 

0 

fa 

fa 

fa 

Eh 

■* 

Eh 

Eh 

* 

Eh 

0 

8 

a 

PA 

0 

a 

0 

o 

O 

ph 

fa 

O 

M 

PM 

o 

0 

0 

0 

+ 

0 

0 

+ 

0 

* 

* 

* 

0 

* 

* 

0 

* 

W 

rH 

rH 

8 

VP 

8 

8 

VP 

a 

X 

(H 

eg 

■VP 

V 

VP 

VP 

SP 

VP 

t> 

> 

> 

li 

H 

1 

It 

a 

1 

i 

Si 

1 

PH 

< 

fa 

W 

fa 

PH 

H 

fa 

X 

EH 

a 

X 

EH 

a 

X 

fa 

a 

< 

fa 

< 

8 

fa 

n 

fa 

fa 

fa 

0 

0 

0 

0 

8 

0 

< 

•fi 

o 

O 

o 

o 

O 

o 

fa 

fa 

fa 

0 

0 

0 

0 

0 

0 

> 

> 

> 

37 


V1“¥2+EK*  ( VBARXI  *COSAXI  -VBARET*COSETA-VBARNU*COSAM5) 


* 

ft 

rH 

^"v 

a 

CM 

CM 

& 

* 

* 

*S‘ 

o 

* 

* 

X 

© 

rH 

* 

rH- 

© 

P 

ft 

ts 

ft 

•> 

•• 

> 

ft 

> 

< 

© 

■< 

a 

< 

Q 

+ 

as 

a 

P 

ft  ft 

rH 

Hr 

a 

Hr 

a 

PH' 

rH 

ft 

ft 

ft 

-ft;  ' 

Eh 

ft 

© 

pH 

© 

ft 

•> 

o 

CM 

o 

Q  ’ 

a 

< 

-ft' 

© 

ft 

•v 

a 

P 

* 

• 

* 

ft  '1 

ft 

* 

a 

© 

a 

a 

X 

ft 

CM 

1 

rH 

-ft  ' 

© 

> 

© 

rH 

eM 

ft 

ph 

«C 

a 

ft 

© 

a 

X 

s 

P 

a 

ft 

© 

h  a 

55 

+ 

© 

a 

i?» 

«>s 

¥%■ 

ffl 

CD 

a 

ft 

© 

O  X 

CO 

a 

O 

’Hr 

a 

ft  ft 

o 

p 

H- 

ft 

HH1 

o 

-•■ 

© 

+ 

Eh  a 

* 

o 

ft 

o 

rH 

&  o 

p 

55 

nr 

* 

> 

< 

-ft' 

* 

e 

ft 

ft 

5S 

-rH 

P 

ft' 

H3 

H  i 

X 

< 

a 

nr 

a 

p 

Eh  ! 

D 

£9 

> 

ft 

© 

Hr 

a  | 

> 

< 

© 

■»* 

ft 

a  j 

53 

+ 

as 

o 

ft 

© 

ft,  i 

p 

nr 

-ft' 

a 

O 

H  i 

H 

Eh 

ft 

* 

a 

ft 

ft 

CQ 

ft 

to 

a 

* 

ft  X 

CO 

a 

ft 

• 

© 

p 

O 

o? 

pH 

© 

© 

i 

H 

o 

a 

© 

i 

© 

o* 

* 

-PTH 

hp 

a 

© 

© 

n  <y 

Eh 

CM 

p-\ 

\ 

© 

Tl« 

a 

co  « 

w 

> 

* 

/PN 

CO 

ft 

a 

a  < 

* 

co 

rH 

© 

0 

© 

CM 

ft 

<3 

ft 

ft 

p 

pH 

a 

X 

£9 

CM 

a 

> 

- 

1 

© 

a 

PH  N 

> 

/H 

CO 

ft 

< 

© 

© 

ft 

i 

a  . 

-H 

+ 

>«< 

• 

M 

* 

H/ 

rH 

ft 

tK 

a 

w 

• 

pH 

a 

• 

* 

* 

© 

rH 

a 

ft  a 

X 

rH 

1 

X 

ft 

ft 

< 

< 

• 

ft 

1 

\  a 

ffl 

i 

CM 

p 

a 

hp 

HP 

nr 

X 

a 

© 

HP 

SR 

00  j 

ft 

CO 

CM 

* 

ft 

p 

H- 

* 

X 

* 

a 

ft 

oo 

rH  j 

B 

o 

* 

* 

* 

ft 

CM 

o 

< 

a 

* 

s 

ft 

t> 

PQ  ( 

a 

a 

* 

< 

* 

* 

* 

CM 

a 

ft 

rH 

tH 

X 

* 

X 

* 

* 

hp 

ft 

ft 

ft 

■> 

ft 

CM 

O  X! 

M 

< 

05 

X 

PN 

r-N 

X 

p 

-*! 

> 

1 

* 

1 

CM 

X 

\ 

ft 

w 

o 

rH 

CM 

•< 

* 

«*! 

V 

< 

X 

W  i 

H 

« 

rH 

rH 

N, 

s 

CO 

* 

* 

ft 

P 

ft 

ft 

X 

ft 

ft 

9V 

*5 

> 

> 

pH 

rH 

t-H 

ft 

* 

ft 

Hr 

ft 

x 

< 

ft 

•3  . 

H 

S-/ 

pH 

w 

OS 

a 

a 

a 

a 

ft 

+ 

X 

PH 

X 

rH 

X 

EH  a 

> 

V 

pH 

ST 

V 

HP 

ft 

a 

PH 

ft 

ft 

ft 

pH 

ft 

ft 

© 

a 

a 

a 

v 

* 

PH 

* 

* 

* 

- — / 

ft 

a 

ft 

* 

a 

* 

X 

ft 

* 

* 

Eh  1 

ft 

* 

• 

X 

• 

pH 

CM 

HP 

ft 

* 

ft 

* 

• 

CM 

ft: 

ft 

«$ 

* 

ft 

& 

B 

X 

CD 

ft 

CO 

rH 

rH1 

ft 

ft 

pH 

li¬ 

CM 

O 

© 

ft 

* 

rn 

S 

ft  a 

w 

X 

* 

X 

PH 

t-H 

E* 

* 

CO 

ft 

• 

© 

© 

ft 

ft 

• 

ft 

PH 

H  ht 

1 

* 

a 

• 

a 

a 

a 

CO 

M 

a 

CM 

© 

• 

ft 

a 

© 

ft 

ft 

P 

CM 

o 

55 

t* 

ft 

ft 

O' 

w 

a 

M 

a 

ft 

o 

It 

a 

© 

a 

I 

O 

> 

II 

M 

II 

a 

ii 

a 

a 

ft 

ft 

ii 

11 

>r 

ft 

ft 

P 

ft 

< 

Eh 

JL 

rH 

ft 

CM 

rH 

CM 

n 

ft 

!i 

* 

ft 

ft 

X 

PH 

it 

ft1 

Hr 

Q 

X 

as 

pH 

II 

pH 

CO 

CO 

a 

« 

pH 

O 

a 

a 

ft 

ft 

ft 

ii 

ft 

ft 

Eh  : 

ft 

pH 

rH 

M 

PH 

PH 

PH 

CO 

CO 

a 

ft 

ft 

> 

a 

ft 

ft 

© 

PH 

PH  | 

r-t 

X 

GO 

X 

X 

a 

a 

CO 

a 

n 

ft 

ft 

< 

a 

a 

ft 

o 

Eh 

QtS  ( 

0 

> 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

a 

a 

P 

p 

ft 

ft 

Q 

ft 

O’ 

P  ft 

a 

r 

■H, 

a 

38 


7HQ  TILDA, 8X, 10HPOS  ST  IMP, 5X, 11HPOS  DYN  IMP/ 4HPSIG, 11X,4HPSIG, 11X, 4HPSIG 
1 IX , 9HMILLISECS,  6X, 4HPSIG, 1IX, 4HPSIG, 11X, 8HPSI-MSEC, 7X, 8HPSI-MSEC/F7 . 2, 


A 


to 

pm 

X 


o> 


. 


to 

Pm 


A 

CSj 


•• 

CO 

Pm 


«S 


ctt 

• 

CO 

Pm 


•4 


cd 

• 

co 

Pm 

a 

X 


CO 

A 

p? 

co 

Pm 

a 

X 

0> 

A 

CS| 


CO 

Pm 


U 


< 

Em 

< 

El 

< 

0 

a 

P 

fa. 

g 

g 

§ 

a 

+ 

+ 

< 

Em 

<c 

X 

o 

05 

pq 

w 

P 

w 

05 

a 

S’ 

Pm 

a 

0> 

H 

p 

a 

o 

0> 

B 

Pm 

0 

+ 

11 

1 

0 

< 

Q 

p 

w 

Em 

“^5 

Em 

+ 

+ 

a 

< 

o 

& 

H 

Pm 

a 

!! 

P 

o 

p 

a 

B 

t 

a 

Pm 

s 

Em 

a 

H 

Pm 

a 

<y 

a 

o 

a 

a 

05 

co 

m 

CO 

c* 

o 

i 


39 


TABLES 
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559.6  152 
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679. C  130 
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6 07. €050 

53.25947 

87*. 130* 
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0 . 3  c 1*70 1 4 
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DATA  FOR  XM103  BLAST  FIELD  COMPUTATIONS 


Parameter 
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PEAK  STATIC  SHOCK  OVERPRESSURE  FOR  XM103 

CANNON  AT  100%  RATED  MAX  PRESSURE  W/O  BLAST  SHIELD 

(Computed  by  assuming  a  point  explosion  at  the 
muzzle  with  50%  energy  utilization.)  qe  -  10°,  h  -  0 
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EXPERIMENTAL  VALUES  OF  PEAK  STATIC  OVERPRESSURE 
FOR  XM103  GANNON  AT  100%  RMP'  W/M0ZZLE  BRAKE  WTV-D8259 


AND  W/0  BLAST  SHIELD  (REFERENCE  12)  QE  -  0° 
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PEAK  STATIC  OVERPRESSURE  FOR  XM103  AT  100% 

RATE©  .MAX  PRESSURE  '»:/  M2A2E2  BRAKE  (V,TV-D8259)  AfJD 
tV/O  BLAST  SHIELD  f  =  0.95,  <*>  =  0,  <?E.=  10°,  h  =  0 
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EXPERIMENTAL  VALUES  OF  PEAK  STATIC  OVERPRESSURE 
FOR  XM1G3  CANNON  AT  100%  RMP ,  W./  MUZZLE  BRAKE 
WTV-F8241  AND  W/0  BLAST  SHIELD  (REFERENCE  12)  QE  -  0° 
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COMPUTED  V A CUES  OF 


PEAK  STATIC. OVERPRESSURE  FOR  XM103  AT  100% 

RATED  MAX  PRESSURE  «/  XM103  BRAKE  (iVTV-F8241)  AND 
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COMPUTED  VALUES*  OF  PEAK  STATIC  OVERPRESSURE  FOR 


XMT03  AT  100%  RATED  MAX  PRESSURE  W/  XM103  BRAKE 


(WTV-F8241)  AND  W/0  BLAST  SHIELD  (*  assuming  point- 
source  decay  from  two  sources)  \>  **  1 . 35,  co  =*  0, 


POSITION  Of  MUZZLE  AT  0°Q£ 
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EXPERIMENTAL  VALUES  OF  PEAK  STATIC  OVERPRESSURE: 

FOR  XMI03  CANNON  AT  100%  RMP  W/MUZZLE  BRAKE 
WTV-F8241  AND  W/0  BLAST  SHIELD  (REFERENCE  12),  QE  =  45° 
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COMPUTED  VALUES  OF 


PEAK  STATIC  OVERPRESSURE  FOR  XM103  AT  100% 

RATED  MAX  PRESSURE  W/XM103  BRAKE  (WTV-F8241) 

AND  W/O  BLAST  SHIELD  *  -  1 . 35,  o> -  0,  QE  -  45°,  h  -  0 
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_* 

COMPUTED  VALUES  OF  PEAK  STATIC  OVERPRESSURE  FOR  XM103  AT  100% 
RATED  MAX  PRESSURE  W/  XM103  BRAKE  (WTV-F8241)  AND  W/0  BLAST  SHIELD 
(*  assuming  point-source  decay  from  two  sources)  b  ■  1.35,  <o  -  0, 
QE  -  45°,  b  -  0 
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POSITION  OF  MUZZLE  AT  0°GLE 


ISOBARJC  PLOT  OF 
COMPUTED  RESULTS  FOR 
XMI03  CAWUOU  AT 
IOO%  RMP  W/XMI03 
BRAKE  AND  W/O  BLAST 
SHIELD  ,  45°  Q.E 
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DISCUSSION  OF  RESULTS 


Comparison  of  the  computed  results  with  the  experi¬ 
mental  results*  presented  in  reference  12,  shows  that  the 
computed  values*  while  displaying  the  proper  directional 
and  decay  characteristics*  remain  consistently  higher 
than  the  experimental  values.  The  discrepancy  in  the  peak 
Static  overpressure  is  nearly  constant  in  magnitude.  This 
fact  allows  one  to  use  the  theory  to  predict  the  "worst ~ 
case’’  or  most  pessimistic  overpressure  that  would  be  en¬ 
countered  in  using  a  particular  muzzle  brake .  The  dif¬ 
ferences  between  theory  and  experiment  may  be  attributed 
to  several  simplifications  in  the  model.  Some  of  these 
are  i 

(1)  The  assumption  was  that  maximum  peak  static 
overpressure  is  developed  at  3"  -»  1  after  discharge  begins. 
This  obviously  depends  to  a  certain  extent  upon  the  abil¬ 
ity  of  the  brake  to  act  as  a  reservoir.  It  was  assumed 
that  conventional  brakes  would  not  appreciably  affect  flow 
rates  into  the  atmosphere  relative  to  those  from  the  bare 
muzzle.  If  this  were  not  so,  one  could  expect  consider¬ 
able  departure  from  theoretical  predictions. 

(2)  Combustion  was  assumed  complete  at  the  muzzle 
with  gases  passing  thru  the  brake  in  a  thermodynamically 
"frozen"  state.  If  this  assumption  were  not  valid* 
departures  therefrom  would  result  in  lower  overpressures 
than  predicted. 

(3)  The  Rankine-Hugonict  relations  used  were  strictly 
one-dimensional,  whereas  the  spherical  shock  was  quasi 
one-dimensional  and  essentially  three  dimensional  in 
character.  Here  again  the  assumption  resulted  in  too 
large  a  predicted  value. 

(4)  Initial  decay  of  the  shock  wave,  i .  e .  at  very 
great  strength,  is  somewhat  more  rapid  than  the  r "3/2 
rate  assumed.  Altho  decay  at  very  low  overpressures 
(below  about  3  psig)  is  not  as  great  as  the  predicted  rate* 
the  total  effect  of  the  latter  assumption  is  to  predict 

a  somewhat  large  value. 

(5)  Constructive  interference  in  the  waves  emanating 

from  the  two  point  centers  of  shock  expansion  was  assumed 
for  all  positions  within  the  shock  field.  This  assumption 
becomes  progressively  less  tenable  as  increases  from 
zero.  In  general,  the  theory  prediets  the  correct  direc¬ 
tional  effects  with  PF~0 ;  however,  it  is  noted  that  the 
theoretical  overpressure  values  for  ■#  23  50°  remain 
somewhat  greater  than  those  at  30°  and  40°  at  all 

corresponding  distances,  D,  from  the  muzzle,  whereas  the 
reverse  situation  is  true  for  the  experimental  values. 
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This  discrepancy  can  no  doubt  be  explained  by  noting 
that  —  contrary  to  assumption  --  complete  constructive 
interference  is  not  likely  for  large  angular  deviations 
from  the  tube  axis. 

In  general,  however,  a  15%  agreement  between  theory 
and  experiment  was  achieved.  Considering  the  complica¬ 
tions,  this  was  all  that  was  expected  and  probably  all 
that  is  required  for  engineering  purposes. 

Since  in  general  best  agreement  with  experiment  is 
achieved  for  PF=0,  this  is  the  value  used  in  the  compu¬ 
tations  for  the  isobaric  plots.  The  fact  that  gaseous 
discharge  is  transonic  rather  than  supersonic  may 
account  for  the  negligible  effect  of  PF.  As  previously 
mentioned,  this  fact  was  suggested  by  the  results  of 
reference  3.  With  PF-O,  this  theory  contains  no 
empirical  constants.  Therefore  the  predictive  ability 
of  this  work  is  no  better  than  the  assumptions,  good 
agreement  being  obtained  when  the  assumptions  are  valid. 

As  can  be  seen  by  inspection  of  the  graphs,  the 
experimental  results  from  the  XM103  with  M2A2E2  brake 
and  with  XM103  brake  bracket  the  computed  results 
obtained  from  the  point  source  explosion  theory. 

This  is  a  fortuitous  result  of  the  choice  Of  brakes . 
Since  the  computed  overpressures  predicted  by  this 
theory  do  not  depend  upon  angular  deviation  from  the 
Z-axis  or  upon  the  momentum  index,  P  ,  of  the  brake, 
quite  severe  inaccuracies  could  result  from  its  use 
as  a  predictive  tool . 

However,  despite  the  differences  in  magnitude 
between  experiment  and  the  results  of  the  point 
source  theory,  the  latter  more  accurately  displays 
the  general  shape  (functional  form)  of  the  decay  of 
static  overpressure  with  distance  than  does  the 
theory  previously  described  in  this  report.  The 
reason  that  this  is  the  case  is  to  be  found  in  the 
assumption  of  a  simple  r“^/2  decay  rate .  Therefore, 
a  modification  of  the  theory  was  made  to  account 
for  a  variable  decay  exponent .  Each  shock  center 
at  maximum  shock  strength  was  considered  to  be  a 
point  source  with  explosive  energy  equal  to  one^ 
half  the  energy  available  in  the  gas  discharged  at 
JT  =  1.  Then,  using  the  pressure  decay  theory 
developed  in  reference  1,  the  pressure  components 
at  the  reference  position  from  the  near  and  far 
shock  spheres  were  computed .  These  were  combined 
as  indicated  in  equations  0.70  and  0.72  to  obtain 
the  measurable  peak  static  overpressure .  Compar¬ 
ison  of  the  computed  results  of  this  modified 
theory  (graphs  7  and  11 )  with  corresponding 
experimental  results  (graphs  5  and  9  )  indicate 


an  extremely  favorable  agreement  and  suggest  the  use  of 
the  modified  theory  as  a  predictor  in  place  of  the 
original  theory. 

It  appears  that  the  single  point  source  explosion 
theory,  badly  underestimates  the  positive  impulse  and 
duration  of  the  positive  phase  and  is,  therefore,  not 
recommended  as  a  predictor  for  these  variables. 
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APPENDIX  I 


The  Rankine-Hugoniot  relations  for  stationary,  normal 
shook  are  given  by  Shapiro  in  reference  18.  Slight  modi¬ 
fications  and  algebraic  rearrangement  of  the  equations 
found  in  Volume  I  on  pages  137  and  following  may  be 
written  here  as  follows.  Notational  conventions  listed 
there  are  preserved  in  what  follows. 

By  requiring  that  Ty  -  Tx  ,  one  obtains 


1.1 


T 

r? 


1  +  -Tx 


H  -  1  »2 


1  + 


-  1 


Also, 


1.2 


T 

Tx 


n  \ 

2 

t  M 

•  y  ' 

“I 

px  j 

i  “x 

1 

l 

From  1.1  and  1.2,  we  have 


1.3 


Ez 

Pt 


X 

*y 


i  +- 


v  -  1 


1/2 


x  +Jl_UL  U2 

1  *  2  7 


From  the  combined  momentum  and  continuity  relations, 
one  can  obtain 

1.4  pJCl  +  M^2)  -  Py(l  +  Ky  My2) 

Elimination  of  p_/p„  between  1.3  and  1.4  yields 

J  X 

M2  (1  +  -  -  m|)  m|(1  +  1  My) 

15  (1  +  **  Mx)2  (1  +  ^y  My) 2 

Calling  the  left  hand  side  of  equation  1.5 


1.6 


F  -  F(MX), 


1-2 


we  have 

.  x  V2 

F  -  (1-2  KyF-  2F) 

i  -  yy  (i  -  2  xy  f)  " 
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Following  Shapiro,  we  let  primed  quantities  refer 
to  conditions  observed  by  a  stationary  observer  in 
region  x  toward  Whom  the  Shock  is  propagating.  Then, 
the  transformation  conditions  relative  to  a  stationary- 
shock  frame  of  reference  are: 


1.19a  a?  -  Z2  a2  £. 

y  Px 


,  where 


b  r 


&x 


ty  Ry 


Using  equations  1.4  and  1.14,  one  can  write 


1.20a  M'  -  x  - 


and 

1.20b  v’  -  a„  M' 

y  y  y 


t  +  y  m2  _  ZZ 

It  ®  X 


vp 


By  definition  1L  *  -JL  ,  or  by  equation  1.17, 

*  a 


1.21 


tg  Yx  Rx  Tx) 


T/5F 


Considering  vx,  the  shock  velocity,  to  be  the 
independent  variable,  one  can  successively  compute  Mx, 

My,  Py,  Ty,  ay.  My,  and  vj  by  sequential 
Px 

application  of  formulas  1.21,  1.6,  1.3,  1.2,  1.19,  1.20a-, 
and  1.20b.  Thus,  tables  can  be  constructed,  functionally 
relating  the  above  variables. 

Utilizing  the  characteristics  of  the  mathematical 
model  of  the  shock  phenomenon,  previously  described,  one 
can  additionally  relate  the  stagnation  temperature  at  the 
center  of  the  shock  sphere,  T°,  to  the  above  parameters, 
By  assuming  that  the  combustion  products  expand  nearly 
adiabatically,  one  can  write 


1 . 22 


£  -  1  +  hjil  Mo 

t2  Z  2 
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2 

However,  it  was  assumed  in  the  model  that  M2  -  1 
and,  therefore,  that  v2  -  a2  -  (g  R2  T2) 1/2  . 

Thus, 

'  *  .  »■  i  i 


.23a  T°  = 

Further, 

(  *2  *  4 

1 . vf  | 

1  2  T 

at  the  ir 

i  *2  g  “2  / 

it  erf ace  v  * 

Finally, 

y 

therefore, 

.23b  T°  -  1 

*2  *  D 

v-2 

.12 

2  X2  & 

R2 

Three  additional  relations  that  were  found  to  be 
useful  to  tabulate  were 


1.24 


vy  =  jv  |Ry  Ty 

V-  V*,  \JWPy 


T° 

1.25  xj_  =  Sr  and 

py 


1.26  ps  =  py  -  p^ 


The  above  tables  were  computed  using  the  following 
constants. 


T 


Poo 


g 

— 

32.17 

»* 

- 

1.4036 

= 

1.3900 

h 

= 

1.2593 

Rx 

= 

53 . 280 

Ry 

= 

53.290 

R2 

- 

70. 036 

Tx 

= 

518.7  °R 

Px 

14.7  psia 
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a, 


(  Vx  g  Rx  Tx) 


v-  Ef  V, 


Rx  % 


Px 


The  variables  X-p  Vy/Vx,  and  ps  are  listed  as 
functions  of  T°  on  pages  40-41  of  this  report . 

For  Vy  =  Yx  =  1.4,  Shapiro  has  derived  the 

relationship  between  the  shock  velocity  and  the  peak 
static  overpressure  shown  below. 

1.27a  <f  -  | 

b  <P  -  Ps/p<*> 


Vx 

2  -  1 

a 

_ 

where 


Also,  for  the  above  conditions,  reference  1. 
the  following  relationship  for  the  dynamic  pressure . 


1.28 


<1 


5  P« 


JL 


2(7  +  <p) 


The  expressions  for  <p  ,  q,  DPP,  SPLSP,  and  DPLSP, 
listed  in  the  EQUATIONS,  were  obtained  by  fitting  data 
given  in  reference  1. 


i 
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APPENDIX  II. 

Derivation  of 
Shock  Model 
Equations 


As  shown  in  figure  2 . la,  flow  of  combustion  gases 
from  a  muzzle  brake  is  conceived  as  issuing  principally 
from  two  ports  of  equal  area.  The  air  which  these 
gases  replace  is  compressed  as  shells  in  front  of  the 
rapidly  expanding  spheres  of  combustion  gases.  The 
spheres  are  shown  here  fully  developed.  We  shall 
return  shortly  to  a  description  of  the  two-shock  model . 

In  figure  2. lb  is  seen  the  single-shock  model. 

The  sphere  of  combustion  gases  is  labeled  region  2; 
the  shell  of  compressed  air  is  labeled  region  y ; 
and  the  surrounding  undisturbed  air  is  labeled  w>  . 

At  maximum  shock  strength*  the  speed  of  the  center 
of  mass  of  the  shock  sphere  in  the  Z-direction  is 
given  by  v2  .  For  momentum  conservation  in  this 
direction* 


2.1 


M 


g 


-  -  jfc 

where  .1  is  the  total  mass  set  in  motion  in  the 
Z-direction.  It  should  be  recognized  that  v„  is  an 
effective  speed  chosen  for  momentum  conservation. 
Considering  only  the  momentum  discharged  from  the 
muzzle,  one  can  obtain  the  actual  average  gas  speed 
from 


g  P(l)  *  P_  v2  -  v? 

M* ( 1 )  Bt  1 . 46881 


Since  the 
occupied  the  volume 
one  can  write 


air  in  region  y  formerly 
7rr|  at  ambient  conditions* 


2.3  a 


.  3 

4/3  TTY: 2 


'X 


4/3  n(K 


r2} 


or 


b  r2  ==  r2(l  v  Vy/Vx) 1/3 
2.4  M* (1)  -  11(1)  + 

TTZ 


Also 
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Now,  the  gas  discharged  to  <?  -  1,  is  considered 
to  occupy  the  volume  4/3  irrx  .  Thus, 


2.5 


2 .6  a 


Trr 


2  _ 


and 


'2  av 


7rr|  -  V2  av 


V°  or 


irr |  «  av  M(1)V°,  where 

_  V, 


av 


2  av 
V° 


2.7 


But,  assuming  ideal  gas  behavior, 
v°  .  *2  T° 


rr 


With  the  additional  assumptions  of  nearly  adiabatic 
flow  within  region  2  and  sonic  flow  at  the  interface. 


2.8 


ta,0  o 

1  *•>.  +  i 

p2  py 

t 

2 

(  V( 


Equations  2.7  and  2.8  yield 

Q  ”2 


Ro  Tq 


2.9a  v 

or 

b  V° 

where 


144  p^  [(  r2  +  l)/2 ] 


(  *2/(  r2  -  D ) 


R2  *1 


144  [ (  T2  +  l)/2  ] 


<  *yc  *2  - 


c  Xi  = 


Finally,  substitution  of  V°  in  2.9b  into  2.6b  gives 


2.10 


47rr? 


av 


R2  *1 


144  (  (  r2  4-  l)/2  ] 


To  obtain  an  expression  for  ^2  av*  we  make  the 
additional  assumption  that  the  radial  velocity  within 
region  2  can  be  written  approximately  as 


2.11  3?-  =  in 

v2  r2 

Thus 

the 

center  t 

>f  mass  frame 

2.12  v  =  v2  J>  , 

wher 

e  0  4 

f<  1  . 

Within  region  ' 

Uadi 

aba tic  i 

flow) , 

2.13  *2  =  ^  - 

1  + 

*2  ”  1 

L . f  \ 

2  J 

\  *2  g  R2  ,T/_ 

V1* 


Also. 


2 . 14a 


*2  g  »2  T 


*2  f 


>  T2  g  R2  Tg  -  V2  g  R9  t°  -|_£2, 

Subtracting  2.14b  from  2.14a  and  using  2.12,  we  have 


2.15 


^2  g  R2  T  -  ^2  g  R2 


+i  h  i 1 W2 


a  *  f) 


By  assumption,  the  mach  number  at  the  interface 
is  unity;  so  that 

2.16  v|  =  ^2  g  R2  t2  Thus, 

2.17  y2  g  R2  T  =  v|  [  1  +  ((  y2  -  l)/2)  (1  -  p)  ]  . 


2.18 


2.13  and  2,17, 
V 


fx  JJ 

'  -  1 
2  A 

f  2 
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l+(( 

!f2-l)/ 

'2)  (1-  p) 

11%1)) 
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By  definition. 


4  it  (1)’ 


2 . 19a  Vay  - 


4  TT  /J(p/v(|>)  )  dj> 

b  Is  (3  /o  |  P2/  1  a)  df 


-1 


Equation  2.19b  was  evaluated  numerically  using  a  value  of 
1.26  for  Kg  .  The  resulting  value  tjg  av  -  1.309  was 

Used  thruout  the  rest  of  the  computations. 

In  the  preceeding  derivations  in  this  appendix,  we 
have  made  use  of  mass  and  momentum  conservation.  Now 
we  shall  appeal  to  the  principle  of  energy  conservation 
to  derive  a  most  important  relationship. 

The  total  energy  available  for  distribution  is  that 
discharged  up  to  S  *  1*  J  H(l) .  Part  of  this  is  involved 
as  translational  kinetic  energy  of  the  gas  relative  to 
the  ground:  m*(1)v|/2  g  . 

Another  part  is  involved  as  stagnation  enthalpy  in 
the  center  of  mass  system: 

2  M(1)t9  .  The  final  major  portion  of  the  energy 

is  invested  in  doing  work  against  the  atmosphere.  This 
work  is  just  the  work  done  in  displacing  the  volume 
4/3  7rr^  against  atmospheric  pressure.  This  portion 


is  thus  144  p0 


4  rr  r2  II  1  +  V 


After  striking  an  energy  balance,  we  have 


2.20 


J  HU)  = 


M*(l)  v|  Yn  R 


2  g 


+  144  P„ 


4  v  r* 


i 1  **i 


Having  written  these  equations  for  the  single  shock 
model,  we  shall  now  examine  what  modifications  must  be 
made  for  a  two  shock  model,  what  must  be  done  to 

account  for  the  presence  of  the  muzzle  brake. 
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The  first  alteration  in  the  above  equations  that 
must  be  made  for  a  model  with  two  shook  spheres  concerns 
the  mass  contained  in  each  of  the  identical  spheres. 
Ignoring  discharge  from  the  front  port*  this  is  one  half 
the  mass  contained  in  the  single  sphere  model.  Thus* 
mass  contained  in  each  sphere  is  M* ( 1 ) / 2 .  Similarly* 
the  momentum  transported  by  each  sphere  is  just  P (l)/2. 
Thus  the  average  speed  of  the  mass  in  each  sphere  is 


g  P(l)/2  -  g  P<1).  as 
H*(l)/2  «*(1) 


before  for  the  single-shock  model.  Thus*  V|  -  Vr>  . 

It  must  be  noted  that  the  total  energy  available  to 
each  sphere  is  now  H(l)/'2.  Now*  by  examining  equations 
2. 10  and  2.20*  it  can  be  seen  that  substitution  of 
M(l)/2*  H(l)/2,  and  vj_  for  M(l),  H(l),  and  respectively 
produces  the  conditions; 


2.21  T°  (two-shock  model)  =  T®  (single-shock  model) 

2.22  radius  of  each  gas  sphere  *; 


4 


By  equation  2.3b  and  2.22* 


2.23 


where  r^  is  the  shock  radius 


for  each  sphere. 

To  predict  the  position  of  each  of  the  fully  developed 
shock  spheres  with  respect  to  the  muzzle  brake*  we  must 
first  apply  momentum  conservation  to  the  system  of  spheres 
and  muzzle  brake.  Regard  figure  2.2  . 

In  figure  2,2a  a  cross  sectional  view  of  a  conven¬ 
tional  two-baffle  brake  is  pictured.  To  simplify  the 
theoretical  treatment*  we  find  the  effective  center  of 
the  side  ports  and  imagine  the  equivalent*  brake  to  be 
pictured  as  in  figure  2.2b  . 

*  A  brake  with  two  side  ports  haying  the  same  axial 
and  vertical  momentum  indices. 


SIDE  PORTS 


For  momentum  conservation  in  the  Z-direction  at  T  «  1, 
2.24  M’1' (1) Vjjj./g  +  Z -impulse  on  brake  at  <!j  -  1  «  Bt(1)  . 

However, 


2.25  Z-impulse  on  brake  (at  ^  -  1)  -  beff  %(!)  . 
Therefore, 


2.26 


'zl 


g  B^xi) 


M 


(1  “  ^ef  f  ^ 


and  from  2.1  and  2.26, 

2.  27a  vzl  -  v2(l  -  *ef  f  >  or 

b  %i  -  n(l  -beff)  . 


Defining 

9.  p.R  X  =  net  Z-axis  impulse  , 

r  .  P  . . "" . . 

Bt 

2.29a  =  -p— (1  -  ^(Esff)  or  ■#  for 

b  Xp  *  1.46881(1  *  Petf)  . 


Thus, 

2.30  v  . 

zl 


V1  Ar 

1.46881 


We  have  defined. 


2.31 


CO  .as 


net  Y-axis  impulse 
P 


1.26  , 


in 
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2.32 


g 


net  Y-axis  impulse  •= 


From  2.2 *.  2.31*  and  2. 32*  one  has 
%  co 

2.33  v  x  =  1.4-6881  * 


Since  is  a  vector  with  components  vxp  vyl* 


and  v2l  , 

_  w  2 

2.34  vxi  »  (vx  -  vyl 


From  figure  2 . 2b*  it  can  be  seen  that 


A  _1  V 

2 . 35a  oC  =  tan  — 


zl 


>  for  vzl  >  0 


b  $  *.  4  ,  f  or  »  1  -  0  i 


c  ^  =  7r  +  tan”1 


*  for  vzl  <  0 


Similarly* 


2.36a  ft  -  tan 


-1  l  ^xl 

V..1 


for  vyl  ^  0  ; 


J-  <  for  Tyl  -  0  • 


Calling  the  distance  from  the  center  of  the  side 
port  of  the  brake  to  the  center  of  one  of  the  shock 
spheres  d-^  *  one  can  write 


2,37  =  net  Z-axis  impulse  and 


3- 


2.38  fl  _  . yl  „ . 

B*f  net  Y-axis  impulse 


by  similarity  of  the  displacement  and  velocity  vectors 
Then,  from  2.28  and  2.37, 


2 . 39 


Also,  from  2  31  and  2.38., 


2 , 40 


d  j  60 

1.46881 


Since  d^  is  a  vector  with  components  x-^  ,  y^^  ,  z^  , 

2.41  xt  -  (df  -  y|  -  zf)1/2  . 

we  shall  now  find  an  expression  for  d-j_  .  Gas  discharge 
Starts  at  the  muzzle.  Therefore,  the  center  of  mass  of 
the  gas  discharged  travels  a  distance  5*  thru  the  brake 
itself  before  reaching  a  side  port.  This  gas  is  expanding 
and  starting  to  Shock  the  Surrounding  air  even  while  in 
the  brake.  Let  us  call  d  the  distance  the  center  of  mass 
of  the  gas  discharged  travels  to  time  of  maximum  shock 
strength,  i.e.  to  -  1  .  Then,  d^  can  be  written  as 

2.42  dx  “  d  -  ... 


The  infinitesimal  mass  discharged  during  the  interval 
dt  is  M  dt  .  The  distance  traveled  by  this  mass  from  the 
time  of  its  ejection,  t  ,  to  the  time  t*  is  given  by 
yav  (t*  ~  t)  ,  where  vav  is  the  average  velocity  of  the 
mass  during  this  interval.  Thus,  the  distance  from  the 
e  center  of  mass  of  all  the  gas  discharged 

Jp*  vav  M  dt  (t*  -  t) 

M(t*) 


muzzle  to  t 
to  t*  is 

2.43  d  - 


Specifically,  we  are  interested  in  t*  such  that 
JT  “  L  Remembering  that  t  =  ^  or  dt  «  d <T  2.43 

CK  —  -  ■"  ■ 

can  be  written  as  1 


2.44 


V  ft  4£  (1  -  S  > 


M(l)  <* 


T 


From  reference  10  we  have 


2.45a  JL.  =  Mt  C4  (C^  +  1) 
T 


Where 


*  «i  -  -*r*  «♦ 


c  ^4  =  y 


y^jr-rj  cc  *  +  i)/2 c  y  -  d) 


d  =  and 


2.46 


0.50688  Mt 


In  our  case  Y  =  =  3- *2®  >  and 


2.47 


d 


^4  ab  init  f1 
0.50688  J0 


(vav/ab  init)  (cl^ 


+  l)"C5(l-J0d^ 


To  complete  our  evaluation  of  d  ,  we  need  an 
expression  for  vav  for  the  interval  (1  -  &  )  .  Let  us 
proceed  by  developing  an  expression  for  the  velocity 
of  the  gas,  Vg  =  vg(  J- )  .  We  note  that  the  velocity 

with  which  the  gas  leaves  the  muzzle  can  be  expressed  as 


2.48 

2.49a 


Vgo  = 

P( 


f-?.C£D  .  From  reference  10  we  have 

(Cr) 

-  Apb  init  f>  >  where 
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b  C 
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r_  =  y  2/i  X  +  IV 

2(y-i) 

3  a  1 

UTk 

2  H 

Also, 


2.50 


A  g 


*  is)  «  - c4 


ab  init 


c  y  +d/2  y 


2. 51 


From  2.48  ,  2.49  ,  and  2.50  , 

Co 
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-gQ- 
b  init 


C4(Cltr  +  1)  ' 


Due  to  atmospheric  entrainment,  the  gas  discharged 
slows  down  as  it  proceeds  away  from  the  muzzle.  After 
traveling  for  a  time  Jr  =  1  ,  we  Shall  assume  that  the  gas 
is  traveling  at  v^  ;  and  that  for  i7"»  1  ,  the  gas 
velocity  asymptotically  approaches  a*.  .  Actually,  of 

course,  for  JT  approaching  <»  ,  vg  approaches  zero,  it  is 
regarded  as  more  realistic  to  assume  that  vg  approaches 
a*,  ,  for  the  purpose  of  establishing  a  function  vg(vT) 
for  the  interval  0  ^  v?"  ^  1  . 

Let 

2.52a  u  =  vg  -  a*„ 


b  u0  “  vgo  ”  *"> 


c  ux  “  vx  -  aoe 


Selecting  a  simple  function  which  passes  thru  the 
points  (vgo,  0)  ,  (Vl,  1)  5  and  approaches  a^  as  ?' 

approaches  <x>  ,  indicates  the  following  form? 

2.53a  u  «  uQ  -  ]_  +  ■*  > 
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b  a  =*  Uq  b 


e  b  ■“  Uq/u-j,  -  1  or 
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d  r  *  s 


q  »  Sj_  *  equation  2.54  becomes 
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By  definition. 
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From  2.56  and  2. 57, 

2-58  «*v  -  -inU  *|j-  ?•> 


Replacement  of  s  ,  etc.  in  2.58  with  their  equival 
ents  (2.52  and  2 , 
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that 

2.60 


Sample  computations  with  typical  Weapons  have  shown 
*1  *  ab  init  * 
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letting  £?'  =  ] 
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1 at ter  expres si on 
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interest,  we 
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7  b  init 


) 


Equations  2.47*  2.51,  and  2 , 61  suffice  to  determine 
oCij  d/ab  init  numerically  as  a  function  of  ab  init/aoo  • 

This  evaluation  has  been  performed  for  ab  iai t/a*,  between 
1.5  and  3. 5,  a  realistic  range  for  the  examples  we  have 


treated.  The  resulting  values  for  cCf  d/ab  init  lie 

between  0.521  and  0.523  .  To  three  significant  figures 
a  value  of  0.523  was  chosen.  Therefore,  equation  2.42 


d  »  0,523  a-b  init  _  ^ 

1  ocT 


2.62 
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APPENDIX  III 


Derivation  of  Shock  Model  Equations: 
Shock  Decay 


As  shown  schematically  in  figure  3.1a,  the  plane 
thru  the  axis  of  the  rear  trunnion  and.  center  of  the 
front  port  of  the  muzzle  brake  has  been  labled  1.  The 
I -plane  is,  in  general,  inclined  at  some  angle  QE  to 
the  H-plane,  a  horizontal  plane  thru  the  trunnion  axis. 

If  the  rear  trunnion  axis  passes  thru  the  tube  centerline, 
the  angle  QE  is  exactly  the  quadrant  elevation;  if  not, 
there  is  some  small,  constant  difference  between  QE  and 
the  quadrant  elevation. 

Using  the  X-Y-Z  coordinate  axes  at  the  center  of 
the  front  port  of  the  muzzle  brake,  the  position  of  the 
center  of  the  near  shock  sphere  is  at  (x, y, z) .  As  seen 
from  figure  3.1,  these  components  can  be  written  as 

3.1a  x  -  Xq+  x^ 

b  y  -  y0+  yx 

C  Z  «  ZQ+  Zj  . 
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In  figure  3.2  we  have  defined  a  reference  position, 
R.  R  is  in  a  plane  normal  to  H  thru  an  axis  parallel  to 
the  trunnion  axis  at  a  distance  f  behind  it .  Within  this 
plane,  R  is  at  a  height  h  above  the  H-plane  and  is  in 
another  plane  normal  to  H,  passing  thru  0,  the  angle  of 
intersection  of  this  plane  with  I  making  an  angle  ■#  with 
respect  to  the  tube  axis.  Note  that  we  have  defined  an 
orthogonal  frame  (  > *  rj  ,  v  )  at  Q  along  ’EPQ.  The 
centers  of  the  near  and  far  shock  spheres  are  labeled 
A  and  6,  respectively.  The  distance  HR  is  called  N  and 
the  distance  BR  is  called  0. 

The  other  given  dimensions  have  been  assigned 
literals  as  shown  in  figure  3.2  *  Also,  from  the  geome¬ 
try  of  the  figure*  we  have  given  some  derived  dimensions. 
In  figure  3.3a  other  derived  dimensions  are  given  in  the 
p ",  h  *  \)  --  coordinate  system  which  serve  to  specify 
the  reference  position  with  respect  to  the  centers  of  the 
near  and  far  shock  spheres . 
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,  and  y-  components 
-  x  I  *  |(G  +  z)  sin|QE) 


With  respect  to  R  the  §  -  > 
of  the  distance  to  A  are:  Pd  sin  v  - 
+  y  cos(QE)  -hi,  and  D  cosi£  +  z  CoStQEj  -  y  sin 
With  respect  to JR  the  fc-  ,  n-  .  and  ,  V-  components  of 
distance  to  B  are:  |  D^simf-  +  x  1  ,  I  (G  +  z)sin(QE)  + 
y  cos(QE)  -  h  1,  and1-  I  D  COS#  +  4  cos(QE)  -  y  sin(QE)  . 

Therefore,  oner  can  write  N  and  0  as  follows.  J 


the 


3.2 


3.3 


N  =  |  (©  sin  #  -  x)2 

( (G  +  z)  sin(QE)  +  y  cos(QE)  -  h)2 
(D  cos  #  +  z  cos(QE)  -  y  Sin (QE) ) 2  J 
0  »  |  (D  sin#'  +  x) 2  + 

+  z)  sin(QE)  +  y  cos 


(D  cos#  +  z  cos  (QE)  -  y  siii(QE))  2  j 


From  figure  3.2,  the  distance,  L,  from  the  center 
Of  the  muzzle  to  the  reference  position  can  be  found  from 


3.4 


L2  =  D2  +  (G  sin(QE)  -  h)2 


Figures  3.3b  and  c  show  the  center~of-mass  velocity 
components  at  A  and  B,  respectively.  The  direction  cosines 
of  the  angles  A ^  ,  A  ^  ,  and  A v  and  their  counterparts 

at  B  are  obtained  by  noting  that  for  the  velocity 
component  projections  along  N  and  0,  the  k-  ,  0-  ,  and 

y -  velocity  components  are  proportional 'to  the  corres- 


ponding  distance 

components.  Thus,  we 

have 

3 . 5a 

cosA^  = 

(D  sin  $ 

'  -  x)/N 

b 

cos  A, |  =* 

(  (G  +  z) 

i  sin(QE)  +  y 

cos (QE)  -  h)/N 

c 

COS A^  = 

(D  cos# 

‘  +  z  cos(QE) 

-  y  sin(QE) )/N 

and 
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3 . 6  a 


b 


cosB 


cos 


? 

*1 


(D  sih$'  +  x)/0 

( (G  +  z)  sin(QE)  +  y  cos(QE)  -  h)/0 


cosB, 


cos 


+  z  cos 


-  y  sin(QE)/0 


To  find  the  center -of -mass  velocity  components  along 
the  ,  rj-  ,  and  y-  axes,  given  the  components  vxl  , 

vyl  ,  and  v2l  *  we  note  that  the  X*  Y,  z  axes  are 

rotated  thru  the  angle  (QE)  about  the  X-axis  to  obtain 
J? ,  fjs  V  axes.  Thus,  one  can  write 


the 


r  .  _  n 

■*  ** 

1  0  0 

?xl 

VSR 

0  cos(QE)  sin(QE) 

Vyl 

0  *  sin(QE)  cos (QE) 

VZ1 

Using  the  direction  cosines,  defined  in  3.5  and  3.6, 
and  the  velocity  components,  defined  in  3.7,  we  find  the 
projections  of  Vj,  in  the  direction  of  the  reference 

position  from  A  and  B  to  be,  respectively, 


cosA^  *•  cosA^  -  v„  cosAjj 


and 


-  v, 


cosB,  -  v„  cosB„  -  Vp  coSBy 


V,  CO&Di  -  V„ 
/  ‘  " 


1 


The  fraction  of  these  velocity  projections  which 
is  effective  in  raising  the  shock  velocity  above  y2  in 

the  ground  fixed  frame  we  have  called  K,  where  K  is 
related  to  the  empirical  propagation  factor,  PF,  by 
the  equation. 

3,8  K  -  PF/1. 46881  . 


Therefore,  one  can  write  the  shock  front  velocities 
at  max  shock  strength  for  the  near  and  far  shock  spheres 
in  the  direction  of  the  reference  position  as 


3.9 

3.10 


vi  **  v2  + 
vf  -  v2  - 


.  cosAt  -  vn  cosAn 

/  F  "  ' 

r,  COSB.  +  Vy.  COSBn 

r  r  '  1 


Vy  cosAy  ) 
Yy  COSBy  ) 


and 
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Using  equation  1,27  with  3,9  and  3.10,  yields 


3.11 

%  “  I  t  <' 

,  v  2 

/|/a)  -  1 

] 

and 

3.12 

?i2  “  i  ( <’ 

i  /  v  2  , 

fj/Si)  -  1 

1  • 

As  mentioned  in  the  discussion  of  the  shock  model,  it 
is  assumed  that  the  overpressure  for  each  shock  sphere 
decays  independently  of  the  other  with  an  r“3/2  dependence 
on  distance.  The  dimensionless  overpressure  components 
at  the  reference  position  due  to  the  waves  from  the  near 
and  far  spheres  are 

I  r.  \  3/2 

3,13  ?31  =  Pll  ( IT  j  and 

I  r-j  V  3/2 

3  ‘ 14  ?32  =  ?12  I  W  ) 

From  the  definition  of  <p  , 

3 . 15  ps$1  *  pm  and 

3.16  ps32  =  ?32  Po® 

The  question  that  is  now  suggested  is:  "How  should 
the  above  overpressure  components  be  combined  so  as  to 
reflect  the  measurable  static  overpressure  at  the 
reference  position?"  The  immediate  answer  for  waves  of 
infinitesimal  amplitude  would  be  the  sum  of  the  two 
components.  For  waves  of  finite  amplitude,  however, 
such  is  not  the  case .  Our  determination  of  the  correct 
functional  relationship  is  as  follows. 


97 


Consider  two  cases.  In  the  first  (figure  3.4a),  a 
single  shook  sphere,  at  maximum  shock  Strength,  is  located 
at  0,  a  distance  L  from  a  reference  position  Q,  The 
radius  of  the  shock  at  maximum  strength  is  R*,  where 
R  <<  L.  We  consider  a  shock  of  such  strength,  <p*,  that 
the  r^3/2  decay  rate  is  tenable.  Then>  the  dimensionless 
overpressure  at  the  reference  position  would  be 

3.17  <p  -*  <p*  (R*/L)3/2  . 

In  the  second  instance  (figure  3.4b),  we  divide  the 
shock  energy  and  volume  at  max  strength,  available  in 
the  first  instance,  into  two  parts,  the  shock  centers 
being  located  symmetricall y  at  a  distance  a  from  the 
line  OQ.  We  choose  a  <<  L.  Under  these  assumptions,  we 
would  not  expect  the  resulting  pressure  at  Q  in  the 
second  instance  to  differ  greatly  from  that  in  the  first . 
Calling  the  equal  overpressure  components  from  the  two 
spheres  at  Q  and  <p2,  one  can  write 

3.18  <p  *  (R/L)3/2  -  f  (  <p2)  ,  where 

f(  ^2)  is  a  function  involving  the  com¬ 
ponents  and  <p2  -  To  preserve  geometric  symmetry, 
we  require 

3.19  f(  fit  f2>  x.  U  f2>  ?1> 

that  is,  the  function  chosen  shall  be  symmetric  with 
respect  to  <p^  and  (p 2  <  Further,  the  function  shall 
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be  non-dimensional ,  One  function  satisfying  these  require 
ments  is 


3,20 


One  might  remark  that,  by  assuming  equal  shock  energy 
and  volume  in  the  first  and  Second  instance,  the  maximum 
peak  Overpressure  in  the  latter  is  identical  to  that  in 
the  former,  as  proven  in  Appendix  II .  Further,  the  .radii 
of  the  spheres  at  max  strength  are  given  by  R*/ (2)1' 3 
Thus, 


3,21  <f1  -  Cj>2  -  M 

V2~ 


Substitution  of  the  latter  values  of  ^  and  2 
into  3. 18  and  3 , 20  gives 


3.22a  f(R/L)3/*  • 


2  MMM 


or 


b  (R*/L)3/2  =  % 


ion, 


;  hence 


3.23a 


-  2 


and 


b  p  -  2  . 

It  must  be  stated  that  this  line  of  reasoning  has 
certain  arbitrary  features  which  from  any  position  other 
than  an  operational  engineering  viewpoint  would  be 
inadmissable.  However,  the  case  treated  in  this  report 
certainly  satisfies  the  assumption  well  enough.  There¬ 
fore,  from  3.20  and  3 . 23b,  one  can  write 


3.24a  cp3  =  (  (pgx  +  <p|2)  and 

b  ps3  "  93  Pop  • 

Finally,  from  equation  1 . 28, 


3.25 


q3 


5Poo  93 

2(7  + 
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